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Abstract 

We generalize the concept of a cycle from graphs to simplicial com- 
plexes. We show that a simplicial cycle is either a sequence of facets 
connected in the shape of a circle, or is a cone over such a structure. We 
show that a simplicial tree is a connected cycle-free simplicial complex, 
and use this characterization to produce an algorithm that checks in poly- 
nomial time whether a simplicial complex is a tree. We also present an 
efficient algorithm for checking whether a simplicial complex is grafted, 
and therefore Cohen-Macaulay. 

1 Introduction 

The main goal of this paper is to demonstrate that it is possible to check, in 
polynomial time, if a monomial ideal is the facet ideal of a simplicial tree. 

Facet ideals were introduced in |Flj (generalizing results in |Vil| and |S VVj 
on edge ideals of graphs) as a method to study square-free monomial ideals. The 
idea is to associate a simplicial complex to a square-free monomial ideal, where 
each facet (maximal face) of the complex is the collection of variables that appear 
in a monomial in the minimal generating set of the ideal (see Definition ^. 4fl . The 
ideal will then be called the "facet ideal" of this simplicial complex. A special 
class of simplicial complexes are called "simplicial trees" f Definition 12. 9|l . The 
definition of a simplicial tree is a generalization of the concept of a graph-tree. 
Facet ideals of trees have many properties; for example, they have normal and 
Cohen-Macaulay Rees rings |F1I . Finding such classes of ideals is in general a 
difficult problem. Simplicial trees also have strong Cohen-Macaulay properties: 
their facet ideals are always sequentially Cohen-Macaulay |F2j . and one can 
determine under precisely what combinatorial conditions on the simplicial tree 
the facet ideal is Cohen-Macaulay |F3| . In |F4| it is shown that the theory is 
not restricted to square-free monomial ideals; via polarization, one can extend 
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many properties of facet ideals to all monomial ideals. All these properties, and 
many others, make simplicial trees useful from an algebraic point of view. 

But how does one determine if a given square-free monomial ideal is the 
facet ideal of a simplicial tree? In Section we give a characterization of trees 
that shows this can be done in polynomial time. This characterization is based 
on a careful study of the structure of cycles in Section|3J The study of simplicial 
cycles is indeed interesting in its own right. In graph theory, the concepts of a 
tree and of a cycle are closely linked to each other: a tree is a connected graph 
that does not contain a cycle, and a cycle is a minimal graph that is not a tree. 
Generalizing to the simplicial case, we use the latter property, together with the 
existing definition of a simplicial tree, to define the concept of a simplicial cycle. 
We then prove the remarkable fact that a simplicial cycle is either a sequence of 
facets connected in the shape of a circle, or a cone over such a structure. This 
in turns yields an alternative characterization of trees, given in Section 

This result enables us to produce a polynomial time algorithm to decide 
whether a given simplicial complex is a tree. The algorithm itself is introduced in 
Section|Sl where the complexity and optimizations are also discussed. Section^ 
focuses on the algebraic properties of facet ideals: in Section 16.11 we discuss 
a method of adding generators to a square-free monomial ideal (or facets to 
the corresponding complex) so that the resulting facet ideal is Cohen-Macaulay. 
This method is called "grafting" a simplicial complex. For simplicial trees, being 
grafted and being Cohen-Macaulay are equivalent conditions |F3j . We then 
introduce an algorithm that checks whether or not a given simplicial complex 
is grafted and discuss its complexity. 

Implementations. The algorithms described in this paper have first been 
coded using CoCoAL, the programming language of the GoCbA system (please 
see http://cocoa.dima.unige.it/). These prototypical implementations can be 
downloaded from |CFS2| . Much more efficient (but less user friendly) C++ 
implementations have been developed for several versions of Algorithm l5.1l using 
the CoCoALib framework (http://cocoa.dima.unige.it/cocoalib/). The C++ code 
is also available at the website |CFS2| . 

2 Simplicial complexes and trees 

We define the basic notions related to facet ideals. More details and examples 
can be found in [FT1IF3] . 

Definition 2.1 (Simplicial complex, facet). A simplicial complex A over a 
finite set of vertices V is a collection of subsets of V, with the property that if 
F E A then all subsets of F are also in A. An element of A is called a face of 
A, and the maximal faces are called facets of A. 

Since we are usually only interested in the facets, rather than all faces, of a 
simplicial complex, it will be convenient to work with the following definition: 
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Definition 2.2 (Facet complex). A facet complex over a finite set of vertices 
V is a set A of subsets of V, such that for all F, G e A, F C G implies F = G. 
Each F e A is called a facet of A. 

Remark 2.3 (Equivalence of simplicial complexes and facet complexes). 

The set of facets of a simplicial complex forms a facet complex. Conversely, the 
set of subsets of the facets of a facet complex is a simplicial complex. This defines 
a one-to-one correspondence between simplicial complexes and facet complexes. 
In this paper, we will work primarily with facet complexes. 

We define facet ideals, giving a one-to-one correspondence between facet 
complexes (or, equivalently, simplicial complexes) and square-free monomial 
ideals. 

Definition 2.4 (Facet ideal of a facet complex, facet complex of an 
ideal). 

• Let A be a facet complex over a vertex set {vi, . . . , v n }. Let k be a field, 
and let R = k[xi, . . . , x n \ be the polynomial ring with indeterminates 
Xi, . . . , x n . The facet ideal of A is defined to be the ideal of R generated 
by all the square-free monomials x i± . . . x is , where {v^ , . . . , Vi s } is a facet 
of A. We denote the facet ideal of A by J 7 (A). 

• Let I = (Mi, . . . , M q ) be an ideal in the polynomial ring k[x\, . . . , x n ], 
where k is a field and Mi, . . . , M q are square-free monomials in x±, . . . , x n 
that form a minimal set of generators for /. The facet complex of I is 
defined to be 5^(1) — {f\ , . . . , F q }, where for each i, Fj = {vj | Xj |Mj, 1 < 
3 < n}. 

From now on, we often use Xi,...,x n to denote both the vertices of A 
and the variables appearing in .F(A). We also sometimes ease the notation by 
denoting facets by their corresponding monomials; for example, we write xyz 
for the facet {x, y, z}. 

We now generalize some notions from graph theory to facet complexes. Note 
that a graph can be regarded as a special kind of facet complex, namely one in 
which each facet has cardinality 2. 

Definition 2.5 (Path, connected facet complex). Let A be a facet com- 
plex. A sequence of facets F\ , . . . , F n is called a path if for alH = 1, . . . , n — 1, 
Fi n F i+ i 7^ 0. We say that two facets F and G are connected in A if there 
exists a path F\, . . . , F n with F\ = F and F n = G. Finally, we say that A is 
connected if every pair of facets is connected. 

Notation 2.6. If F, G and H are facets of A, H G means that H n F C 
GnF. The relation <f defines a preorder (reflexive and transitive relation) on 
the facet set of A. 

Definition 2.7 (Leaf, joint). Let F be a facet of a facet complex A. Then F 
is called a leaf of A if either F is the only facet of A, or else there exists some 
G e A \ {F} such that for all if e A \ {F}, we have H ^ F G. The facet G 
above is called a joint of the leaf F if F n G ^ 0. 



3 



It follows immediately from the definition that every leaf F contains at least 
one free vertex, i.e., a vertex that belongs to no other facet. 

Example 2.8. In the facet complex A = {xyz, yzu, uv}, xyz and uv are leaves, 
but yzu is not a leaf. Similarly, in A' = {xyu,xyz,xzv}, the only leaves are 
xyu and xzv. 




Definition 2.9 (Forest, tree). A facet complex A is & forest if every nonempty 
subset of A has a leaf. A connected forest is called a tree (or sometimes a 
simplicial tree to distinguish it from a tree in the graph-theoretic sense). 

It is clear that any facet complex of cardinality one or two is a forest. When 
A is a graph, the notion of a simplicial tree coincides with that of a graph- 
theoretic tree. 

Example 2.10. The facet complexes in Example 12. 81 are trees. The facet com- 
plex pictured below has three leaves F\, F2 and F3; however, it is not a tree, 
because if one removes the facet F4, the remaining facet complex has no leaf. 




The following property is proved in |F3I Lemma 4.1]: 

Lemma 2.11 (A tree has two leaves). Every tree with two or more facets 
has at least two leaves. □ 

3 Cycles 

In this section, we define a simplicial cycle as a minimal complex without leaf. 
This in turns characterizes a simplicial tree as a connected cycle-free facet com- 
plex. We further show that cycles possess a particularly simple structure: each 
cycle is either equivalent to a "circle" of facets with disjoint intersections, or to 
a cone over such a circle. 

Definition 3.1 (Cycle). A nonempty facet complex A is called a cycle if A 
has no leaf but every nonempty proper subset of A has a leaf. 
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Equivalently, A is a cycle if A is not a forest, but every proper subset of A 
is a forest. If A is a graph, Definition 13.11 coincides with the graph-theoretic 
definition of a cycle. The next two remarks are immediate consequences of the 
definitions of cycle and forest: 

Remark 3.2. A cycle is connected. 

Remark 3.3. A facet complex is a forest if and only if it does not contain a 
cycle. 

In the remainder of this section, we provide a complete characterization of 
the structure of cycles. 

Definition 3.4 (Strong neighbor). Let A be a facet complex and F, G 6 A. 
We say that F and G are strong neighbors, written F ~a G, if F ^ G and for 
all H e A, F n G C H implies H = F or H = G. 

The relation ^a is symmetric, i.e., F G if and only if G ~a F. Note 
that if A has more than two facets, then F ^ A G implies that F n G ^ 0. 

Example 3.5. For the facet complex A' in Example 12.81 xyu ^a> xzv, as 
their intersection x lies in the facet xyz. However, xyz ~a' xzv and similarly 
xyz ~a' xyu. 

Remark 3.6. Suppose A is a facet complex, and A' C A. Let F,G £ A'. If 
F ~a G, then F ~a' G. The converse is not in general true. 

Remark 3.7. We have F ^a G if and only if G is strictly maximal with respect 
to on A \ {F}, i.e., for all H ^ F, G H implies G = H. This is a simple 
restatement of the definition. 

It turns out that a cycle can be described as a sequence of strong neighbors. 
The following lemma follows directly from Definition 13.41 

Lemma 3.8. If A is a facet complex with distinct facets F,G\,G2 such that 
F ~a G\ and F ~a G%, then F is not a leaf of A. 

Proof. If F is a leaf, there exists a facet H ^ F such that G\ H and 
G2 H, which by Remark l3.7l implies that G\ — G2 = H, a contradiction. □ 

Corollary 3.9. Let IS. be a facet complex, and let F%, . . . , F n be distinct facets 
with n ^ 3, such that F\ ^a F2 ~A • ■ ■ ^A F n ~a ^i- Then . . . , F n } has 
no leaf. 

Proof. This follows directly from Remark 13.61 and Lemma l3~%l □ 

Lemma 3.10. Suppose A is a facet complex and F,G £ A. If F is a leaf of 
A \ {G}, but not a leaf of A, then F ~ A G. 

Proof. Suppose H is some facet such that F nG C H, but H ^ F and H 7^ G. 
Since F is a leaf for A \ {G}, there exists a facet H' € A \ {G} such that 
L n F C H 1 for all L € D \ {F, G}, and so F (1 H C ff'. But now we have 
FnGCFnHCH', which implies that F is a leaf of A, a contradiction. □ 
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Proposition 3.11 (A cycle is a sequence of strong neighbors). Suppose 
A is a cycle, and letn = \A\. Then n 3, and the facets of A can be enumerated 
as A — {F\, . . . , F n } in such a way that 

Fx ~a F2 ~a ■ ■ • F n ~a -Fi, 

and Fi ^ a Fj in all other cases, so that each facet is a strong neighbor of 
precisely two other facets. 

Proof. First note that since A is not a forest, n ^ 3. We begin by showing that 
each facet has at least two distinct strong neighbors. Let F E A be a facet. 
Since A is a cycle, A \ {F} is a tree. The subset A \ {F} also has cardinality at 
least two, and therefore has two distinct leaves, say G and FL, by Lemma \2. Ill 
Since neither G nor H are leaves of A (because A is a cycle), we have F ~a G 
and F ~a H by Lemma T3.1UI 

Now we can simply choose F\ arbitrarily, then choose F2 ^ F\ such that 
Fi ~a F2, then for every i ^ 3 choose Fi such that F$_i ~a Fi and Fi ^ 
Fi-i, Fi-2- Since A is finite, there will be some smallest i such that Fi = Fj for 
some j < i. Then A' = {Fj, . . . , Fi-i} has no leaf by Corollary 13.91 so A' = A. 
It follows that j = 1 and i — 1 = n. Finally, suppose that Fk ~a Fi for some 
k ^ / — 2, where k > 1 or I < n. Then {Ft, . . . , Fk, Fi, . . . , F n } has no leaf by 
Corollary 13. 91 contradicting the fact that it is a tree. □ 

The converse of Proposition 13 . 1 II is not true. 

Example 3.12. The facet complex A is not a cycle, as its proper subset A' 
(which is indeed a cycle) has no leaf. However, we have Fi ~a F2 ~a G 
F3 ~a Fi ~a Fi, and these are the only pairs of strong neighbors in A. 





F2 








F2 




Fl 




F3 


A' = 


Fl 




F3 




F4 








F4 





Lemma 3.13. If A is a cycle, written as F± ^a F2 ■ ■ ■ ~A F n ~a Fi, then 
for each i, Aj = A \ {Fi} is a tree with exactly two leaves Fj_i and Fi+x, with 
joints Fi-2 and Fi+2, respectively. 

Proof. We know that Aj is a tree, so it has at least two leaves. By Lemma f3. 81 
Fi—i and Fi+i are the only choices. By Remark |3.7I Fj-2 is the only possible 
joint for Fi—i, and is the only possible joint for Fi+i. □ 

The following lemma will be fundamental for the classification of cycles. 

Lemma 3.14. Let A be a cycle with facets f / G e A. If F ^a G, then 
F n G C H for all H € A. 
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Proof. We first prove the claim in the special case where F ~a H. Indeed, since 
F is a strong neighbor of exactly two facets, there must be some L ^ G,H such 
that L F H. Then Lemma 13.131 implies that H is a joint of F in the 
tree A \ {L}, and therefore F n G C if, or equivalently, F j^g 

Now consider the general case. By Proposition 13. Ill the facets of A can 
be enumerated as F\ ~a F 2 ~a ■ ■ • ^A F n ~a F\. Assume, without loss of 
generality, that F — F% and G = Fi, where 2 < i < n. By repeated applications 
of the special case above, we have 

F <g F 2 • ■ ■ Fi-l- 
In the other direction, we similarly have 

F < G F n ^ G F n -! ■ ■ ■ F +1 . 

Therefore, F (1 G C Fj for j = 1, . . . , n. □ 
Lemma 3.15. Let A be a facet complex, and let 

A = P| F and A' = {F \ A | F £ A}. 

FeA 

Then A' is a facet complex. Moreover, A is a cycle if and only if A' is a cycle. 

Proof. For each F € A, let F 1 = F \ A. Since A is a facet complex, we have 
F % G for any two distinct facets F, G £ A, which clearly implies F' % G' . So 
A' is a facet complex. Let F be any subset of A, and let V = {F' | F £ A} be 
the corresponding subset of A'. Then for any triple of facets F,G,H £ F, we 
have F G F' ^h' G' . Therefore, T has a leaf if and only if V has a 

leaf. It follows that A is a cycle if and only if A' is a cycle. □ 

Theorem 3.16 (Structure of a cycle). Let A be a facet complex. Then 
A is a cycle if and only if A can be written as a sequence of strong neighbors 
F\ ^a i*2 ■ ■ ■ ~A F n ~a F\ such that n 3, and for all i,j 

n 

Fi n Fj = P) F k if j ^ £ - + 1 (mod n). 
fe=i 

Proof. Let A be a cycle. Then by Proposition 13 . 1 II and Lemma f3. 141 A can be 

written as a sequence of strong neighbors with the desired properties. 

Conversely, suppose that A is written as a sequence of strong neighbors 
Fx ~ A F 2 ~a • • • ~a F n -a Fi such that F, n Fj = flfc=i F k V j ? i - 
1,£,£ + 1 (mod n). By Lemma 13.151 we can without loss of generality assume 
that HLi F k = 0- 

By Corollarv l3.9l A has no leaf. Suppose A' is any nonempty proper subset 
of A. We need to show that A' has a leaf. Suppose F{ £ A' and -Fi+i ^ A'. 
There are two cases: 
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1. Fi-i i A'. In this case, since Fi C\ F k = $ for all F k e A' \ {F^, Fi is a 
leaf. 



2. Fi_i e A'. In this case, F, n F fc C for all F k G A' \ {FJ, and so 
is again a leaf, this time with Fi-% as a joint. 

So A is a cycle and we are done. □ 

The implication of Theorem l3.1tjl is that a simplicial cycle has a very intuitive 
structure: it is either a sequence of facets joined together to form a circle in 
such a way that all intersections are pairwise disjoint (this is the case where the 
intersection of all the facets is the empty set in Theorem 13. 16(1 . or it is a cone 
over such a structure fLemma l3.15l) . 

Example 3.17. The facet complex A is a cycle. The facet complex T is a cycle 
and is also a cone over the cycle V . 




4 Characterization of trees 

We now consider the problem of deciding whether or not a given facet complex 
is a tree. We refer to this problem as the decision problem for simplicial trees. 

Note that the naive algorithm (namely, checking whether every non-empty 
subset has a leaf) is extremely inefficient: for a facet complex of n facets, there 
are 2 n — 1 subsets to check. Also note that the definition of a tree is not inductive 
in any obvious way: for instance, attaching a single leaf to a tree need not yield a 
tree, as Example l2 . 1 01 shows . This seems to rule out an easy recursive algorithm. 

Nevertheless, we demonstrate that the decision problem for simplicial trees 
can be solved efficiently. This is done via a characterization of trees given in 
this section. 

Definition 4.1 (Paths and connectedness outside V). Let A be a facet 
complex, and let V be a set of vertices. We say that a sequence of facets 
Hi, ... , H n € A is a path outside V in A if for alH = 1, . . . ,n — 1, (Hi <lHi + i)\ 
V ^ 0. We say that two facets F, G € A are connected outside V in A if there 
exists a path Hi, ... , H n outside V in A such that Hi = F and H n = G. 

Note that in case V = 0, this coincides with the definition of connectedness 
from Definition 12. 51 



8 



Notation 4.2. If F, G%, G 2 are three distinct facets of A, then we define A^ 1 ' 2 
to be the following subset of A: 

a g^g 2 = { H e A \ HnF = G 1 nG 2 }U {Gi, Gal- 
Definition 4.3 (Triple condition). Let A be a facet complex. We say a triple 
of facets (F,Gx,G2) satisfies the triple condition if G\ ^f G 2 and G2 Gi, 

G G 

and if G% and G 2 are connected outside F in the facet complex A F X ' 2 . 

We note that the definitions of A F 1,G2 and the triple condition have changed 
from an earlier version of this article |CFS1| : they have been simplified. 

Example 4.4. Consider A in Example 13. 121 Then the triple (FbF^Fj) sat- 
isfies the triple condition. This is because F4 F2 and F2 ^Fi F4. More- 
over Ap 2 ' Fi = {F2, F3, F4, G}, and a path connecting F2 and F4 outside F\ is 
F 2 , F 3 , F 4 . 

However, (G,F 2 ,F 3 ) does not satisfy the triple condition, since F 2 ^g F3 
(and F 3 F 2 ). Also A^ 2 :F3 = {F 2 ,F 3 }, and F 2 and F 3 are not connected 
outside G. 

Proposition 4.5 (A triple is part of a cycle). Let A be a facet complex. A 
triple (F, G±, G2) satisfies the triple condition if and only if there exists a cycle 
A' C A such that F, G\, G 2 <E A' and G x ~ A , F ~ A / G 2 . 

Proof. Suppose (F, Gi,G 2 ) satisfies the triple condition. Then by definition, 
G\ ^f G2 and G2 ^f Gi. Choose a minimal (with respect to inclusion) path 
Hi, ... , H n outside F that connects Hi = Gi to H n = G2 ■ Note that minimality 
implies that for j > i + l, {H t C\H.j)\F = 0. We claim that A' = {F, Hi, ... , H n } 
is a cycle with 

F ~ A , Hi ~ A , . . . H„ ~ A , F. (1) 

(a) F ~A' Gi and F ~ A / G 2 . 

G G 

If F n Gi Q Hi for some i, 1 < i < n, then since -f/^ € A^ 1 ' 2 , we have 
F n Gi C Hi n F = Gi n G 2 C G 2 . This implies that Gi ^ F G 2 , a 
contradiction. So F ~ A / Gi, and similarly F ~ A / G 2 

(b) Hj ~ A , for i = 1, . . . , n - 1. 

Since {H t n flt+i) \ F 7^ 0, we have that If, n Hj+i % F. By minimality 
of the path, if Hi D -ffi+i C iJj for some j > i + l, then iJj D Hj+i C 
i/i n i?j C F, a contradiction. The case j < i is similar. 

This shows Q. To finish the proof that A' is a cycle, we must show that it 
meets the remaining condition of Theorem 13. 161 If n = 2, there is nothing to 
show; assume therefore that n ^ 3. By definition of Ap U ° 2 , F n H j = Gi n G 2 
for j = 2, ... ,71 — 1, and so PIgsA' G = Gi n G 2 . Also, if j > i + 1, then 
Hi C\ Hj C F by minimality of the path, therefore 

ff< n Hj = n F) n (f^ n F) = G t n G 2 = n GeA , G. 
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So A' is a cycle. 

Conversely, suppose that A' is a cycle containing F, G\ and G 2 , written as 
F —A' G\ —A' Hi —A' ■ • • ~A' H n —A' G 2 ~A' F, where n ^ 0. 

From the strong neighbor relations it follows that G\ ^f G 2 and G2 G\. 
It also follows that the above sequence of strong neighbors provides a path from 
Gi to G2 outside F. We only need to show that for i = 1, . . . ,n, Hi C\ F — 

d n G 2 . 

If A' = {F, Gi, G2} we are done. So assume that n ^ 1. 

We know i?j f\ and so by Lemma EH H t n F C d n G 2 . On 

the other hand, since G\ ^A' G2, Lemma E>.14l implies the opposite inclusion 
Hi n F 2 Gi n G 2 . It therefore follows that Hi H F = Gi n G 2 and we are 
done. □ 

An immediate implication of Proposition ^. 51 is an (algorithmically) efficient 
criterion to determine whether or not a facet complex is a tree. 

Theorem 4.6 (Main Theorem). Let A be a connected facet complex. Then 
A is a tree if and only if no triple of facets in A satisfies the triple condition. 



5 A polynomial- time tree decision algorithm 

By Theorem 14.61 to check if a facet complex A = {Gi,...,G;} is a tree, we 
only need to check the triple condition for all triples of elements of A. The 
checks themselves are straightforward. Since the triple condition for (F, G, G') 
is clearly unchanged if one switches G and G' , we can limit triple checking to the 
elements of the set {(F, G t , Gj) e A 3 | G, ^ F ^ Gj,i < j}. The procedures 
for the basic steps follow immediately from the earlier definitions. 

Algorithm 5.1 (Tree decision algorithm). 

Input: a connected facet complex A = {Gi, . . . , G/} with n vertices. 
Output: True if A is a tree, False otherwise. 

1. For each triple (F, G, G') € {(F, G h Gj) E A 3 | G t ^ F ^ G 3 , i < j} 

(a) If G G' or G' ^p G, continue with the next triple. 

(b) Build A%' G ' . 

(c) If G and G' are connected outside F in A F ' , return False. 

2. Return True. 

The correctness of this algorithm is an immediate consequence of Theo- 
rem The algorithm uses very little memory; the input A requires nl bits, 
and A F ^ C A requires I bits. The memory required to perform the connected- 
ness check and to store the various counters is negligible. Thus, memory locality 
is good, and the computations can generally take place in the cache. 
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Remark 5.2. In the process of checking the triple condition for a triple (F, G, G 1 ) 
that is part of a cycle, we build a connection path outside F. Clearly, any such 
path can be reduced to a minimal connection path {-Hi, . . . , H n } outside F for 
G, G', and therefore, by the proof of Proposition 14.51 {F, H%, . . . , H n } forms a 
cycle. Therefore, an easy modification of Algorithm 15 . II allow us to produce the 
set of all the facets F £ A that are part of some cycle, and a cycle A' F 3 {F} 
for each of them. 

5.1 Complexity 

For each triple it is trivial to see that steps (a) and (b) can be performed with 
cost 0(n) and 0(nl) respectively. For step (c), the following holds. 

Lemma 5.3. Let A be a facet complex with I facets over n variables such that 
F,G,G' are distinct facets of A. The connectedness outside FofG,G'£A can 
be determined with time cost 0(nl). 

Proof. First of all we substitute A with the set {H\F \ H £ A}. We then define 
n + 1 equivalence relations Po, . . . , P„ on the set {1, . . . , I}. Po is the identity 
relation, i.e., each equivalence class is a singleton. For each j = 1, . . . , n, consider 
the vertex Vj and the set Xj — {i \ Vj £ Fi}. Let Pj be the smallest equivalence 
relation such that Pj-i C Pj and such that for all i, i' £ Xj, £ Pj. Then 

facets Fi and F^ are connected if and only if (i,i') £ P n . With a suitable data 
structure for representing equivalence relations, the complexity of the procedure 
above is 0(nl). □ 

Consequently, step (c) of the tree decision algorithm can be performed at 
cost 0(nl). Thus, the total complexity of the tree decision algorithm is as 
follows: in the worst case we have to check 3 • (3) = l Ski£i±k-£l = 0(l 3 ) triples. 
The complexity of the steps (a)-(c) is 0{nl) and hence the total complexity of 
the algorithm is 0(nZ 4 ). 

Example 5.4. Consider the facet complex A = {xy, xz, yz, yu, zt}. We have 
to check 3 • (3) = 30 triples. We start with the triple (xy, xz, yz). 

• xz ^ xy yz since xy n xz — x % y — xy n yz. Similarly yz ^ xy xz. 

• xz and yz are connected outside xy in the complex A%y ,yz = {zt, xz, yz}. 

We have hence discovered that A is not a tree. A more unlucky choice of 
facets could have brought about the checking of 27 useless triples before the 
discovery that A is not a tree, the other two useful triples being (yz, xy, xz) and 
(xz,xy,yz). 

Example 5.5. Some statistics for a bigger random example. Consider the facet 
complex A = {Ika, qik, tykj, wuv, rjb, eioab, gdc, zv, rtj, qrvm, gzm, tgzb, 
rgvm, qlav, qeocn, ikfaz, bn, ekjs, pfvn, wtodv}. We discover that it is not 
a tree after checking 4 facets; we performed the connectedness check only once. 
If one checks all 3 • ( 2 3 °) = 3420 triples, one finds that 445 of them require a 
connectedness check, and 403 of them reveal that A is not a tree. 
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Example 5.6. The facet complex {2^2^+12^+2 | i = 1, . . . ,400} is trivially a 
tree. Checking this by a direct application of Algorithm 15.11 requires dealing 
with 3 • ( 4 °°) = 31,760,400 triples, and takes about 12.6 seconds on an Athlon 
2600+ machine for our C++ implementation. All the timings in the remainder 
of this paper refer to this machine. 

5.2 Optimization 

The runtime of Algorithm 15.11 can be improved by introducing some optimiza- 
tions. First, note that if F is a facet such that no triple (F, G, G') satisfies the 
triple condition, then by Proposition 14.51 F cannot be part of any cycle of A. 
Therefore, F can be removed from A, reducing the number of subsequent triple 
checks. We refer to this optimization as the removal of useless facets. 

Example 5.7. We check the tree {2^2^+12^+2 | i = 1, . . . , 400} of Example 15. 61 
with a version of Algorithm 15.11 with removal of useless facets. This requires 
checking 10, 586, 800 triples and takes about 3.46 seconds. 

An important special case of a "useless facet" is a reducible leaf, as captured 
in the following definition: 

Definition 5.8 (Reducible leaf). A facet F of a facet complex A is called a 
reducible leaf if for all G, G' G A, either G ^ F G' or G' ^ F G. 

A reducible leaf is called a "good leaf" by Zheng Z . 

Remark 5.9. The facet F is a reducible leaf of A if and only if F is a leaf of 
every A' C A with F G A'. 

The remark immediately implies that a reducible leaf cannot be part of a cy- 
cle. Thus, it can be removed from A, and the algorithm can then be recursively 
applied to A \ {F}. We were not able to find a tree without a reducible leaf; 
in fact, Zheng |2j conjectured that this is always the case. Checking whether 
a given facet F is a reducible leaf requires ordering all facets with respect to 
^f, which takes 0(nllogl) steps. A reducible leaf can thus be found in time 
0(nl 2 log/). Therefore, if Zheng's conjecture is true, the tree problem can be 
decided in time 0(nl 3 log/). But even if the conjecture is not true, remov- 
ing all reducible leaves at the beginning of Algorithm 15. II is still a worthwhile 
optimization. 

5.3 Optimization for sparse complexes 

Let A be a facet complex with / facets. If every F <G A intersects a substantial 
(ps /) number of facets, then the number of cycles is probably high and our 
algorithm is usually able to detect one of them easily. If this does not happen, 
we can exploit the "sparseness" of the facet complex in our algorithm. 

For the remainder of this subsection, A will be a facet complex with / facets 
over n vertices such that the maximum number of neighbors of a facet F G A 
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is d and the maximum number of vertices of a facet F € A is v. Note that trees 
are the hard cases for our algorithm, since all the triples have to be checked. 
Also note that, if A is a tree, then / ^ n. This follows by induction on I, from 
the fact that every leaf contains at least one free vertex. 

5.3.1 Connection set algorithm 

To check if A is a tree it is sufficient to check the connected triples only. For 
each facet F (I facets): first construct the set of all facets G connected to F 
(called the connection set, at cost 0(lv)), then for all G,G' in the set (O(eP) 
pairs) perform the triple check on (F,G,G') (cost 0(nl) per triple). We call 
this optimization of Algorithm 15 . II the connection set algorithm. The total cost 
is 0(nl 2 d 2 y The space required to construct the connection sets is O(d), hence 
negligible. If the complex is not sparse (d sa I, v n), the complexity is 
the same as Algorithm 15. II However, for sparse examples, this optimization is 
clearly worthwhile: 

Example 5.10. We check the tree {2^2^+12^+2 | i = 1, • • • , 400} of Example l5.6l 
with the algorithm detailed above. We deal with 398 triples and spend 0.2 
seconds. 

Example 5.11. The facet complex {xiXi+i ■ ■ ■ 2^+200 | i = 1, • • • , 3200} is a tree 
but not sparse. Tree checking with the connection set algorithm is still quite ef- 
ficient; it requires dealing with 61, 013, 400 triples, and takes about 140 seconds. 
Without any optimization, the number of triples to check is 16, 368, 643, 200 and 
the time spent by the algorithm is > 2 days. 

5.3.2 Incidence matrix algorithm 

The connectedness relation for a facet complex A can be represented by a graph 
through an incidence matrix. This matrix can be built and used during the 
tree checking algorithm. Since creating incidence matrices from a complex is a 
relatively expensive operation, we build them in steps, exploiting at every step 
the relations already computed. 

We compute the connectedness relation for A at cost 0(l 2 d). Then for every 
facet Fe Awe compute the "connectedness outside F" relation for A, at cost 
0(nld). Then for every triple (F, G, G') (there are 0(<i 2 ) of them) we compute 

the "connectedness outside F" relation for A G ' G at cost 0(dv + Id). Using 

this additional structure, we do not actually need to build A G,G , and we can 

check connectedness outside F in A G,G using the connectedness relations at 
cost O(ld). We call this optimization of Algorithm 15.11 the incidence matrix 
algorithm. 

The total complexity for this algorithm is hence 0(nl 2 d + ld 3 v + l 2 d 3 ^j . If A 
is not sparse (v n, d I), then this algorithm has roughly the same complexity 
as Algorithm 15. II 

On the other hand, if d « v i=a \fl ~ \fn,, which is a reasonable assump- 
tion for sparseness, then the complexity of the incidence matrix algorithm is 
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while the complexity of the connection set algorithm is 0(/ 4 ) and 
that of Algorithm O is 0(l 5 ). 

6 Algebraic properties of facet ideals 

We now study facet ideals from a more algebraic point of view. In particular, we 
are interested in ways to determine whether a given facet complex A is Cohen- 
Macaulay, meaning whether R/J 7 (A) is a Cohen-Macaulay ring. We first need 
to introduce some new terminology. 

Definition 6.1 (Vertex covering number, unmixed facet complex). Let 

A be a facet complex. A vertex cover for A is a set A of vertices of A, such 
that A n F for every facet F. The smallest cardinality of a vertex cover of 
A is called the vertex covering number of A and is denoted by a(A). A vertex 
cover A is minimal if no proper subset of A is a vertex cover. A facet complex 
A is unmixed if all of its minimal vertex covers have the same cardinality. 

Example 6.2. Consider the two facet complexes in Example 12.81 We have 
a(A) — 2. Also, A is unmixed as its minimal vertex covers {x,u}, {y,it}, 
{y,v}, {z,u} and {z,v} all have cardinality equal to two. We further have 
a(A') — 1, but A' is not unmixed, because {x} and {y, z} are minimal vertex 
covers of different cardinalities. 

The following observations are basic but useful. 

Proposition 6.3 (Cohen-Macaulay facet complexes FIT IF3 | ~> ■ Let A be 

a facet complex with vertices in x± , . . . , x n , and consider its facet ideal I = J^(A) 
in the polynomial ring R = k[x\, . . . , x n \. Then the following hold: 

(a) height / = a(A) and dim R/I = n — a(A). 

(b) An ideal p = [x^ , . . . , Xi a ) of R is a minimal prime of I if and only if the 
Set {Xi^ , . . . , Xi J is a minimal vertex cover for A. 

(c) Ifk[x%,... ,x n ]/ T(A) is Cohen-Macaulay, then A is unmixed. 
6.1 Grafting 

One of the most basic ways to build a Cohen-Macaulay facet complex is via 
grafting. 

Definition 6.4 (Grafting [F3 ). A facet complex A is a grafting of the facet 
complex A' = {G\, . . . , G s } with the facets Fi, . . . , F r (or we say that A is 
grafted) if 

A = {F 1 ,... !J F r }U{G 1) ... ! G s } 
with the following properties: 

(i) Gi U . . . U G s C Fi U . . . U F r ; 
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Figure 1: Three different ways of grafting the facet complex A. 

(ii) Fi, . . . ,F r are all the leaves of A; 

(iii) {G 1 ,...,G s }n{F 1 ,...,F r } = 0; 

(iv) For i f j, Fi n Fj = 0; 

(v) If d is a joint of A, then A \ {Gi} is also grafted. 

Note that the definition is recursive, since graftedness of A is defined in terms 
of graftedness of A \ {Gi}. Also note that a facet complex that consists of only 
one facet or several pairwise disjoint facets is grafted, as it can be considered 
as a grafting of the empty facet complex. It is easy to check that conditions (i) 
to (v) above are satisfied in this case. It is also clear that the union of two or 
more grafted facet complexes is itself grafted. 

Example 6.5. There may be more than one way to graft a given facet complex. 
For example, some possible ways of grafting {Gi, G2} are shown in Figure ^ 

The interest in grafted facet complexes, from an algebraic point of view, lies 
in the following facts. 

Theorem 6.6 (Grafted facet complexes are Cohen-Macaulay [F3J). Let 

A be a grafted facet complex. Then ^(A) is Cohen-Macaulay. 

Even more holds when A is a tree. 

Theorem 6.7 ((EH Corollaries 7.8, 8.3]). If A is a simplicial tree, then the 
following are equivalent: 

(i) A is unmixed; 
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(ii) A is grafted; 
(Hi) !F(A) is Cohen- Macaulay. 

6.2 Graftedness algorithm 

A direct application of Definition l6.4l is not very convenient for checking whether 
a given facet complex A is grafted, since at each step of the recursion, one 
potentially needs to check condition (v) for several of the Gi, and this leads 
to a worst-case exponential algorithm. In order to arrive at a more efficient 
algorithm, we characterize graftedness as follows: 

Lemma 6.8 (cf. | F3l Remarks 7.2, 7.3]). A facet complex A is grafted if 
and only if (1) for each vertex v, there exists a unique leaf F such that v G F. 
and (2) all leaves of A are reducible. 

Sketch of the proof. First, assume that A is grafted. Condition (1) follows from 
(i), (ii) and (iv). The fact that all leaves are reducible is shown by induction on 
the number of facets of A. The converse is also shown by induction. Suppose A 
satisfies (1) and (2), and let {Ft, . . . , F r } and {Gt, ■ • • , G s } be the sets of leaves 
and non-leaves, respectively. Conditions (i)-(iv) hold trivially. Further, if Gi 
is a joint, then Ft, . . . , F r are still reducible leaves of A \ {Gi} by Remark 15.91 
Also, there are no additional leaves in A \ {Gi}, since none of the Gj have 
free vertices by Condition (1). Therefore, A \ {Gi} satisfies (1) and (2) and is 
therefore grafted by induction hypothesis, proving (v). □ 

The algorithm for checking if a facet complex is grafted follows immediately 
from Lemma 16.81 

Algorithm 6.9 (Graftedness algorithm). 

Input: A facet complex A with I facets and n vertices. 
Output: True if A is grafted, False otherwise. 

1. Build the lists T = {Ft, ...,F k } (leaves of A) and Q = {Gt, G m } 
(facets of A which are not leaves). 

2 - H Ucee G £ er F > return False. 

3. If 3 F , F' E T such that F n F' ^ 0, return False. 

4. If 3 F e T that is not a reducible leaf, return False. 

5. Return True. 

6.3 Complexity 

The leaf checking cost is 0(nl), hence the cost of step 1 is 0(nZ 2 ). The cost 
of steps 2 and 3 is 0(nl). For step 4, there are k facets F to check. Check- 
ing whether F is reducible takes 0(nl log I) steps as mentioned in Section l5~2l 
Therefore the total cost for step 4 is 0(nl 2 logl), and this is the cost of the 
algorithm. 
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Example 6.10. Let A = {xyz,yzu,ztu,uv,tw}, with T = {xyz,uv,tw} and 

Q = {yzu, ztu}. Then UggS ^ — Ufgjp ^ = i 2 "' 2/> z ' u ' w ' w } anc ^ 

a;yz n to = w n tit; = 0. Additionally, we check that each F e JF is a reducible 

leaf by showing that the set {F n G \ G S 5} is a totally ordered set under 

inclusion. For example, if F = xyz, then this set is equal to {yz, z} which 

is totally ordered. This holds for all F E J 7 , and hence the facet complex is 

grafted. 
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